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Abstract 

A quantitative theory is developed, which accounts for heating artifacts in 
three-pulse photon echo (3PE) experiments. The heat diffusion equation is 
solved and the average value of the temperature in the focal volume of the laser 
is determined as a function of the 3PE waiting time. This temperature is used 
in the framework of nonequilibrium spectral diffusion theory to calculate the 
effective homogeneous linewidth of an ensemble of probe molecules embedded 
in an amorphous host. The theory fits recently observed plateaus and bumps 
without introducing a gap in the distribution function of flip rates of the 
two- level systems or any other major modification of the standard tunneling 
model. 

I. INTRODUCTION 



Glasses are commonly considered as nonequilibrated, structurally disordered solids. The 
amorphous network performs configurational fluctuations on a wide range of time scales. At 
low temperatures, tunneling through conformational barriers is assumed to be the dominant 



relaxation mechanism. This view is equivalent to the two-level system (TLS) mo that 
has successfully accounted for many of the anomalous low-temperature properties of glasses.ll 
In the standard formulation, it is assumed that TLS have on a logarithmic time scale a broad 
distribution of energy splittings, E, and relaxation rates, R, 

P(E, R) dEdlogR ~ constant x dEdlogR. (1) 

As a consequence, the observed relaxation dynamics in glasses depends on the experimental 
time scale. This has far-reaching implications for the linewidth of an optical transition of 
a probe molecule embedded in a glassy host: its transition frequency becomes a dynamical 
quantity.! This phenomenon is commonly referred to as spectral diffusion (SD). It prevents 
a definition of a homogeneous linewidth as in crystalline materials. 

In an optical experiment, the absorption spectrum of an ensemble of chromophores doped 
into an amorphous host is measured. Due to inhomogeneities in the local strain or electric 
fields experienced by individual chromophores, this spectrum (the so-called inhomogeneous 
line) is several orders of magnitude broader than the linewidth of a single molecule. Hence, 
to monitor the local structural glass dynamics, line-narrowing techniques like two-, three- 
pulse (stimulated) photon echoes (2(3)PE) and hole burning (HB) have to be used.i~i Being 
conceptionally similar, the main difference of these techniques is the sensitivity to relaxation 
processes on different time scales: fast dephasing processes (pico-nanoseconds) caused by 
rapid TLS flips and vibrations of individual chromophores! are measured in 2PE; linewidths 
measured in 3PE (nano-milliseconds) and HB (microseconds-several days) are additionally 
broadened due to SD. The rates of SD may vary over many orders of magnitude. By 
convention, the linewidth measured in a 2PE experiment is called homogeneous. In the 
standard model, Eq. ([I]), a linear increase with temperature is predicted. The linewidth 
measured in 3PE and HB experiments is commonly referred to as "effective" homogeneous 
linewidth. According to Eq. (|1|), this line also broadens linearly with temperature. By the 
same token, SD induces a logarithmic broadening as a function of the waiting time (time 
between the second and third pulse in a 3PE, or between burning and reading a hole in a 



HB experiment 

Very recently, there has been broad interest in 3PE of Wiersma's group© ) - '- 6 ), which — 
contrary to the standard model — showed a plateau between 10~ 6 and 10~ 3 s. They in- 
vestigated the chromophore zincporphine in deuterated ethanol glass (ZnP/EtOD). Meijers 
and Wiersma interpreted their data by postulating a gap in the distribution of TLS rates 
in the microsecond regionJiiN~( e ) Zimdars and Fayer emphasized, shortly after that, the 
contradiction between hole burning data of Fayer's group@ and the 3PE data of Wiersma 
group. El( 6 ) _ ( e ) To explain this difference, they invented a model, which is based on the 
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sumptions that (1) the coupling between the chromophores and the glass can vary as a 
function of the TLS relaxation rates, and (2) TLS with relaxation rates in the microsecond 
domain couple to the chromophores only via an electric dipolar interaction. According to 
this, non-polar chromophores such as those used in Wiersma's experiment will not sense a 
perturbation by those electric TLS resulting in a gap in the relaxation rate distribution. 
The relevance of this model has been discussed by Silbey and co- worker s&® The same 
authors could fit Meijers' data qualitatively without assuming a gap, but starting out from 
a modified distribution function, which predicts a slower-than-logarithmic line broadening 
in the microsecond region .@( a )~( fe ) Silbey's model is based on computer simulations of a 
Lennard- Jones model glass.0 Finally, it was Zilker and Haarer who shed light into the dis- 
cussion by showing in the first 3PE experiment below 1 K that the line broadening does not 
only level off in the microsecond region, but in fact decreases again after several hundred 
microseconds.0 They investigated the chromophore zinc meso-tetraphenylporphine in poly- 
methylmethacrylate (ZnTPP/PMMA), a polymer glass. As the authors pointed out, this 
reversibility indicates laser induced sample heating artifacts, which may also be responsible 
for the gap previously observed in Wiersma's group. 

It is the content of this work to clarify this point by treating the laser induced heating 
effect in a quantitative way, and reanalyzing the 3PE data of Meijers and Zilker. 

The paper is organized as follows: in Sec. II, we brievely describe the experimental 



apparatus used in Ref. |Tg and give an estimation of the heating effects for EtOD and PMMA 



using the details of Meijers' and Zilker's experiment; in Sec. Ill, we present a simple model 
for the time dependence of the temperature in the focal volume of the laser beam by solving 
the heat diffusion equation for appropriate boundary conditions; in Sec. IV, we calculate 
the optical linewidth in the framework of nonequilibrium spectral diffusion^; in Sec. V, we 
compare with Meijers' and Zilker's 3PE data, and, in Sec. VI, we conclude. 

II. EXPERIMENTAL 

A. Stimulated photon echoes 

The stimulated photon echo measurements in Ref. [18| (see Fig. 1) were performed with a 
tunable dye laser (Coherent 702-1CD) synchronously pumped by a mode-locked argon-ion 
laser (Coherent Innova 200-10). The output pulses (6 ps width) are cavity-dumped and 
amplified in two dye cell amplifiers which use a two-stage design. Each chain is pumped by 
its own frequency-doubled nanosecond Nd:YAG laser (Spectron SL404G) at a repetition rate 
of 10 Hz. After each stage a saturable absorber is used to suppress amplified spontaneous 
emission. Two beams with equal intensities are obtained from the resulting 1 /iJ pulse of the 
first amplifier by a beam splitter. Beam 2 is fixed in time, whereas beam 1 can be delayed 
via a motorized delay line (DL1). The third pulse can be electronically delayed by picking 
any pulse from the 76 MHz pulse train emitted from the dye laser. Synchronization between 
the dye and the YAG lasers is accomplished via the reference signal of the mode-locker 
electronics, which triggers a digital delay generator (Stanford Research Systems DG535). 
The latter controls the firing of the flash lamps and the Q-switches of the YAG lasers and 
is also used to trigger the cavity-dumper of the dye laser. This design allows us to create 
pulses which can be electronically delayed with respect to each other from 100 ns up to 
100 ms. An optical delay line (DL3) is used to realize waiting times shorter than 7 ns. The 
optical density of the sample was about 1 at the wavelength of the study, 597 nm. The laser 
was focused on the sample to a spot size of approximately 100 /iin; the pulse energy was 
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about 3 nJ. The good optical quality of the samples enabled us to perform measurements 
of stimulated photon echoes without using an optical upconversion scheme.ll The signal 
is spatially filtered to block most of the scattered light, and separated from the excitation 
beams by diaphragms. Finally, it is detected by a photomultiplier tube and registered with 
a boxcar integrator. The echo intensity depends nonlinearly on the energies of the excitation 
pulses. This makes an A/B-measurement for compensation of excitation energy fluctuations 
impossible. Therefore only data points, for which the energy of each laser pulse — measured 
by the two photodiodes (D) — is within a ±10 % window, are processed by the computer. 
Several scans of the delay line were averaged for each signal curve. 

Zinc meso-tetraphenylporphine (ZnTPP, Porphyrin Products Inc.) was dissolved at a 
concentration of 1.5 x 10~ 4 M in the monomer methylmethacrylate; the polymerization of 
the mixture was initiated at 50° C. The resulting samples had a thickness of 1 mm and a 
radius of 4 mm. The fluorescence lifetime 7\ of ZnTPP in PMMA was measured as 7.5 ns 
in a transient-grating experiment .@ 

The sample was mounted in a 3 He bath cryostat (Janis Research Co.) providing a 
temperature range from 0.4 up to 4.0 K. The temperature was measured with a calibrated 
Germanium resistor (Lake Shore GR-200A-100) attached to the copper sample holder. Tem- 
perature control was achieved with an accuracy of ±0.01 K by heating the internal charcoal 
pump of the cryostat. 

B. Heating effects in PMMA and EtOD 

Let us assume that at t = the total sample is in equilibrium with the helium bath 
at a temperature T and illuminated by the (first and the second) laser beam. The laser 
irradiation excites via So — > Si transitions those chromophores of the sample which are in 
its focal volume Vf OC ai = 7ra 2 L p . The latter is a cylinder given by the spot size a ~ 50 /iin and 
the penetration depth L p « 0.5 mm of the laser beam. The sample is cylindrically shaped 
with a radius Rq of 4 mm and thickness L = 1 mm. The chromophores are assumed to be 
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homogeneously distributed in the illuminated volume, and the laser beam propagation is 
taken to be perpendicular to the circular surface of the sample. The excited chromophores 
irradiate heat, Q, over the fluorescence lifetime t/i ~ 1 — 10 ns (in PMMA: Tfi ~ 7.5 ns, 
in EtOD: r/j ~ 2.7 ns) due to intersystem crossings, S\ — > T, to the triplet state. Since 
metalporphines such as ZnP and ZnTPP undergo these transitions with probabilities of 
more than 90 % ISC yield, up to 25 % of the incident laser energy can be transformed into 
heat during this process. If the heat release is considerably faster than the heat diffusion 
process out of the focal volume, we can assume that after a time Tfi the chromophores have 
heated up the entire focal volume to a maximum temperature T\ which is determined by 
the relation 

Q= [ Tl c(T)dT. (2) 

In PMMA, the specific heat is given by§ c(T) = 4.6T/iJ/gK 2 + 29.2 T 3 /iJ/gK 4 . Un- 
der the conditions of Zilker's 3PE experiment^ (pulse energy of 3 nJ, ISC yield of 
97 %, A£(singulet - triplet) = 3800 cm" 1 ) for ZnTPP/PMMA, Q = 20 //J/g is released 
into the sample which results in a local temperature increase inside the focal volume of 
AT = T\ — T ~ 0.55 K for T = 0.75 K. At T = 1.5 K the heating effects amounts into 
AT ~ 0.2 K, and no heating effect is expected at 3 K. 

In EtOD glass, the specific heat is not known below 2 K. The value at two 2 K is roughly 
by a factor two larger than in PMMA: CEtoo(2 K) ~ 350 /iJ/gK.0 The chromophore used in 
Meijers' experiment, ZnP, has very similar optical properties as ZnTPP, which was used in 
Zilker's experiment: ISC yield of 98 %, A.E(singulet — triplet) = 3700 cm _1 .0 The applied 
pulse intensity in the 3PE experiment of Meijers and Wiersma has been quoted to be "less 
than 200 nJ per pulse" in Ref. 15(a) and "less than 100 nJ per pulse" in Ref. |15|(d). Our 



estimation of the heating temperature, T 1; is based on a pulse intensity of 100 nJ/pulse. 
Since the optical properties of both chromophores are nearly identical, one can estimate 
the heat release into the sample in Meijers' experiment by multiplying Zilker's result by 
100/3, which amounts in Q ~ 650 /uJ/g. In PMMA, this provides a temperature increase of 



AT ~ 1.6, 0.8, and 0.6 K at To = 1.75, 2.4, and 3 K, respectively. Since the specific heat in 
EtOD is roughly by a factor 2 larger than in PMMA, one has to divide these numbers by 2, 
yielding T x ss 2.55, 2.8, and 3.3 K at T ss 1.75, 2.4, and 3 K, respectively. 



III. HEAT DIFFUSION 

Our goal is to determine the heat flow out of the focal volume by solving the heat diffusion 
equation 

^T(r,t)- J DV 2 T(r,t) = (3) 
for the local temperature T(r,t), where the diffusion constant 

D = - (4) 

pc 

is controlled by the heat conductivity k, the mass density p, and the specific heat c 
of the sample. To calculate the line broadening, the spatially averaged temperature 
T(t) = (T(r, t)) rg y focal will be, eventually, plugged in the formalism of nonequilibrium spec- 
tral diffusion. The underlying assumption is that the phonons are in quasi-thermal equilib- 
rium with the momentary temperature T(t) in the focal volume, i.e., the TLS relaxation 
rates R{T) are given by R[T(t)}. This is reasonable, because the phonon relaxation time is 
much shorter than the diffusion times (|^) — (0). The replacement T(r,t) — > T(t) is clearly 
approximate. In a correct treatment, the r dependence of the temperature has to be handled 
on the same footing as the r dependence of the dipole-dipole interaction (~ r~ 3 ) between 
the chromophore and the TLS in the host, which would, however, render the model too 
complicated. 

There are three time scales for the heat diffusion process: the radial and longitudinal 
diffusion time over the time scales 

r a = a 2 jD, (5) 
r L = L 2 p /D, (6) 
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and the time scale over which heat radiates into the helium bath [note c(He) 3> 
c(PMMA/EtOD)] 



which is controlled by the Kapitza surface resistance between a solid and liquid helium 



Although k and c depend on temperature, we will assume in the following that D is tem- 
perature independent. Orders of magnitude for r a , t l , r r are given for PMMA in Table 
I. 

Stimulated echo experiments are performed on the time scale between 0.1-1 ns to 10- 
100 ms. The plateau or the bump appears in the data between about 10 /is and 100 /is. This 
suggests the following picture: For t < r r ~ 2 ms heat cannot radiate into the helium bath 
due to the large value of the boundary resistance R K . Accordingly, the probe itself operates 
as a heat bath and absorbs the heat flow out of the focal volume in direction radial to the 
laser beam. At tl i~ a any inhomogeneous illumination of the focal volume in longitudinal 
direction becomes unimportant. The radial flow sets the time scale for equilibration with 
the helium bath after about 10 — 100 /is. At this time scale, we would expect a bump in 
the stimulated echo data. After several hundred microseconds, we further expect a return 
to equilibrium spectral diffusion at T = T . Triplet heating can be excluded, since r r is 
much less than the time period between two consecutive echo pulse series (~ 75 ms). Both 
experiments confirm this picture (cf. Fig. 2 and 3). 

Let us now perform our analysis more quantitatively. To solve the diffusion equation, 
appropriate boundary conditions have to be posed. First, to keep the model simple we 
neglect any heat conduction in propagation of the beam. This might be a crude approxima- 
tion, however being justified by r a <C r r and because the laser beam penetrates almost the 
whole sample (penetration length L p ~ 0.5 mm). The heat conduction in direction radial 
to the beam is governed by diffusion into the sample. Hence, we have to solve the diffusion 



r r = RkPcV{ ocsl \, 



(7) 



Rk — 



0.05 TK 4 m 2 



(8) 



7ra 2 T 3 W 
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problem in an infinite (boundary effects due to a finite i?o = 4 mm are not important in this 
experiment) cylindrical medium in which heat is produced by point sources (chromophores) 
at a rate ~ AQ(r)e _t//rfi (r = \/x 2 + y 2 ). The amount of heat, AQ(r), released by the chro- 
mophore at position r is governed by the laser intensity, I(r), at that point in the sample. 
We assume that the intensity profile across the diameter of the laser beam is Gaussian, 

-r 2 /2a 2 

where P is the total amount of beam power (in W) used. The exact solution to this problem 
is very complicated and beyond the scope of this paper. Instead we will treat a heat diffusion 
problem, which is equivalent in all physical aspects to the problem mentioned above. We 
assume that the chromophores have heated the focal volume to a higher temperature 7\ 
after the time t^. If <C r a , the temperature profile is then well approximated by the 
laser intensity profile, 

AT(r,0) = 2ATe~ r2/2a2 . (10) 

The normalization, we have chosen, is such that the average temperature in the focal volume 
at t = is AT, i.e., AT = J °° dr 2irrAT(r, 0) p- r2 l 2a2 /2ttcl 2 . The temperature at a later time, 
t, is given by integration over the heat kernel^ 



1 /-oo /-oo p-r 2 /2[a(t)] 2 

AT(r, t) = -±- dx' dy' AT(x', y', 0) e -K— ') 2 Hv~v'?]/^t = 2 AT (11 
AirDtJ-00 J-00 o" t ar 



where 



a(t) = V2Di + a 2 . (12) 

Hence the average temperature AT(t) = / °° dr 2tiv AT(r, t) p~ r2 / 2a2 /2iia 2 at a later time, t, 
reads 

mt) = T^fj-- (13) 

1 + t/T a 

Eventually, to account for the delay of the heat release over the time scale Tfi, we subtract 
e - t / T fi from Eq. ([131) . This corresponds to assuming that heat is produced at the rate 

9 



(Qp/ T fi) e i,/r/i P er un it time and unit volume yielding a final temperature increased AT = 
Q/c. This, finally, gives 

T(t) = To + (Ti - To) {y^- " e-'^j . (14) 

Of course, this equation is only an approximation for the real solution of the diffusion problem 
posed above. As a consequence, the numerical value for AT and r a obtained from the fits by 
using fll^D have to be seen as an order of magnitude estimation. However, our model contains 
all physical characteristics: (i) heat is supplied over the time scale t/j with an exponential 
rate, (ii) the temperature returns algebraically (non-exponentially) to equilibrium after the 
time scale r a . We have checked that (i) and (ii) are independent of the initial temperature 
profile (pT0|) . For instance AT(r, 0) oc 0(a — r) gives the same kind of behavior. 

IV. OPTICAL LINEWIDTH 

Hu and WalkerE and Suarez and SilbeyEH have shown that three-pulse photon echo ampli- 
tude decays as exp[— (Fi(r) +T 2 (r, £))tls] due to TLS-flipping in the amorphous host. Here 
(F 1 (t))tls is the two-pulse echo decay (or dephasing decay), and (T 2 (r, £))tls is the decay 
that depends on the separation, t, between the second and the third pulse. (...)tls denotes 
the ensemble average over the TLS-relaxation rates, R, and the TLS-energy splittings, E. 
If r lies well within a hyperbolic distribution of relaxation rates R, i.e., -R m i n < 1/r < -R ma x 
and P(R) ~ 1/-R, the effective, i.e., t depending dephasing rate [1/T 2 ](i) is defined by 
(Fi(t) + F 2 (t, £))tls — 2r/T 2 (t). Analogously, the two-pulse echo decay rate, l/T 2i2PE is 
defined by (Ti(t))tls = 2t/T 2i2 pe. Subtracting the lifetime contribution defines the pure 
dephasing rate 1/T 2 * 2PE = 1/T 2i2 pe — l/2Ti. The effective homogeneous linewidth measured 
in a 3PE experiment is defined by 

[l/7rT*](t) = [1/ttT 2 *] 2PE + [1/ttT 2 *] sd (0 (15) 

where the last term, arising from T 2 (r, t), is the contribution due to spectral diffusion. For 
TLS and phonons in thermal equilibrium at temperature To with a hyperbolic distribution 
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in TLS-relaxation rates, R, and a flat distribution in the TLS-energy splitting, E, (cf. Eq. 
([[])) one finds the well-known result 

[l/7rT*](t) =KT [3M + \og(R eS t)} + r PLM . (16) 

Here, K is a collection of constants proportional to the dipole-dipole interaction strength 

between the TLS and the chromophore, and R e s is an effective relaxation rate averaged 

over the TLS splittings, E, (see below). The second term, [l/7rT 2 *] SD (/f:) oc log(R e fft), is the 

waiting time dependent contribution due to spectral diffusion. The first and the last term 

describe pure dephasing. At very low temperatures the term linear in T is dominant. This 

term arises from Fi(t) and is characteristic for TLS-induced dephasing in amorphous solids. 

The last term, characterized by an activated temperature dependence 

exp(-AE/k B T) 
" b [l-exp(-AE/fc B T)]2' (17) 

also appears in molecular mixed crystals. It is a pure dephasing mechanism arising from vi- 
brational modes of individual chromophores (pseudo-local phonon modes).! This mechanism 
usually starts to be dominant in the temperature regime between 1 and 3 K. 

Recently, Silbey et alS^ have introduced a new model based on molecular dynamics 
simulations on a Lennard- Jones computer glass.0 In the simulations for Ni-P at a given 
composition, it was found that the distribution function P(A) of the tunneling frequency, 
A, is hyperbolic for A/ks less than ~ 1CT 3 K, but as A increases, P(A) is best fitted by 
a form 1/A 1_l/ where v increases to ~ 0.2 as A/k B increases beyond 1 K. With a typical 
maximum relaxation rate of about 10 10 s _1 at 1 K, this implies that the standard model 
[y = 0) does not apply for waiting times shorter than ~ 10~ 3 — 10~ 4 s. Although these 
results are not quantitatively accurate for all glasses, we will see that they provide a more 
reliable picture in the investigation of heating effects in the 3PE data below. The distribution 
function of TLS-relaxation rates and energy splittings is modified in this model according 
to 

P{E, R) = constant x gWjWT! (lg) 
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(We set Silbey's second parameter fi = 0& a ^) This yields a two-pulse decay varying as 
exp[— (2t/T 2> 2pe) 1 ^^ 2 ] ; which is very close to an exponential decay for fits of experimental 
data. In the limit v <C 1, the effective homogeneous linewidth is then found to behave like 

\i/'nT;}{i) = KX +v (Qv + (2M[i-(Rest)- u/2 }) 

+ r PLM , (19) 

where 

R -u/2 = I E — dE E v sech 2 (E/2k B T) R m ^ 
cS J^dEE»sech. 2 (E/2k B T) ' 1 ' 

and (replacing E max /2k B T by infinity) 



K V = K dxx v sech.\x) (21) 



o 



with 0^=0 = 3.66 and 0^=0.25 ~ 3.4 (note K U=Q = K). 

To calculate the spectral diffusion contribution in the case of heating over the time scale 
t, it is more convenient to use the formalism of Black and Halperin.0 According to this, the 
contribution to the optical linewidth due to spectral diffusion is determined by the number 
of TLS, fif(t), which have flipped an odd number of times during the time t. Denoting the 
initial occupation of the upper (lower) state by n + (0) (ri_(0)), and the upwards (downwards) 
relaxation rate by W + {W_), a simple master equation approach yields 

n-WMW-^^ (22) 
R 

where R = W + +W^. In our case, all TLS are initially in thermal equilibrium with the helium 
heat bath at T . Hence, n + (0) = [exp(E/k B T ) + l] -1 and n_(0) = [1 + exp(— E /k B To)Y l . 
If we assume quasi-thermal equilibrium for the phonons, i.e., if we assume that the phonons 
are at the momentary temperature (|b4]), the relaxation rates become time dependent and 
satisfy W+(t)/R(t) = [exp(E/k B T(t)) + l]" 1 and W_(t)/R(t) = [1 + ex V (-E /k B T{t))}- 1 , 
with the one-phonon rate 

R[r,T{t),x] = rR msx [T(t),x], (23) 
i2max[T(t), x] = cT V coth[xT /T(t)], (24) 
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where x = EjlksTQ and c ~ 10 10 K~ 3 s -1 . In the experimentally relevant parameter regime, 
the exact numerical solution of the master equation with time dependent rates is very exactly 
described by the following analytical expression: 

n/(t) = (1/2){1 - tanh[£/2£; B To] tanh[E /k B T(t)}} (l - e -«W0,*n . (25) 



The average of over the TLS-energy splitting, E, and the dimensionless relaxation 

rate, r = i?/i? max , is now performed using the distribution function (|T8|) . With this, the 
effective homogeneous linewidth at the momentary temperature T(t) reads 

[l/7rT*](t) = K u T^ +u [& u + f u (t,R max [T(t)})} + r PLM , (26) 

where 



/ °° dx [1 - (t JU TO,!])^ [l-taph( ar) tanh(xT /T(t))] 
(i//2)/ °°dx ^sech 2 (x) 



/i/(t, -rtmaxi-t (tJJJ = 77~77^ roo 3 ; 2/_.n • \ 2 '> J 



The result for the standard model arises in the limit v — > 0. In Fig. 2 and 3, [l/7rT 2 *](t) 
[(|19D and (H)] is shown together with 3PE data of Ref. P(d) for ZnP in EtOD, and Ref. @ 



for ZnTPP in PMMA. In Fig. 4, 2PE data of ZnTPP and tetra-tert-butyl-terrylene (TBT) 
in PMMA are shown together with fits made with 



l/vrT 2 *] 2PE = K v T^B, + r PLM . (28) 



V. COMPARISON TO EXPERIMENT 

Meijers and WiersmaEl have studied the 3PE decay of various chromophores in a number 
of different glasses. In order to fit their data they were forced to modify the standard model 
(|I|) by postulating the existence of a gap in the distribution of TLS flip rates on the time 
scale between micro- and milliseconds. For the system ZnP/EtOD, they assumed that the 
decay rate is linear in log(t) from the earliest times (~ 1CT 10 s) to 1CT 6 s, flat until 1CT 3 s, 
and then once more linear above 10~ 3 s, with the same slope as in the early time regime. 
As discussed in Ref. |l6|(a), the physical nature of this gap remains totally unclear. 
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Zilker and HaareJll have studied for the first time a 3PE decay below 1 K. They found 
at 0.75 K for the system ZnTPP/PMMA that the decay rate increases from the earliest 
times (~ 10~ 10 s) to 10~ 6 s, levels off, and decreases after ~ 300 /is. As the authors noted, 
this reversibility clearly indicates laser heating artifacts in the sample. Clearly, the question 
arises whether laser heating might also turn out to be responsible for the observation of a 
plateau in Meijers' and Wiersma's data. 

In Ref. ^(a)-(b), Meijers' and Wiersma's data have been analyzed with the modified 
standard model using Eq. ([19]). The authors have found qualitative agreement with better 
quality fits emerging for the high temperature data. Zilker and HaarerEl also analyzed 
successfully their data at 1.5 and 3 K, but not at 0.75 K, within this model. 

In Fig. 2(a), we show the 3PE data for ZnP/EtOD at three temperatures, taken from 
the thesis of Meijers^ 1 -^ with two fits using Eq. (|T9|), dashed dotted line, and Eq. (p6|), 



full line. As noted already in Ref. the need for a gap in the distribution function has 
disappeared. As shown in the figure, the predicted value of the 2PE decay rates also agree 
with experiment. One clearly sees from this plot, that including heating artifacts into the 
theoretical description improves the fits to the 1.75 and 2.4 K data, significantly. We have 
demonstrated this more clearly in Fig. 2(b). The initial temperature increase due to laser 
heating is AT = 0.9, 0.4, and 0.3 K at T = 1.75, 2.4, and 3 K, respectively. Both 3PE and 
2PE show consistently that the contribution to pure dephasing of pseudo-localized librational 
modes of the guest chromophores in the amorphous matrix is insignificant even at 3 K for 
ZnP in EtOD glass.S^ 

In Fie;. 3(a ), the 3PE data for ZnTPP/PMMA taken from Zilker and Haarer@ are plotted 
along with our theory in the same way. Again the data can be fitted with no gaps in the 
distribution function, and with parameter values close to those used in Figs. 2(a)-(b). The 
reversibility of the line broadening for the 0.75 K data is well reproduced by our theory, as 
depicted in Fig. 3(b). This clearly shows the influence of laser heating artifacts. No heating 
effect is observed at 3 K data; the temperature increase at 0.75 and 1.5 K is found to be 0.75 
and 0.2 K, respectively. However, the predicted value of the 2PE decay rates is much lower 
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than the experimental value. The 3PE data for ZnTPP/PMMA can well be fitted without 
(or with only a small, cf. Fig. 4) contribution from localized phonon modes. Contrary to 
this, the 2PE data clearly show a deviation from a power law temperature dependence above 
~ 2 K, which indicates active localized vibrations of the chromophores. In Fig. 4, we have 
illustrated this using recently published 2PE data for ZnTPP/PMMA@ together with new 
results for TBT/PMMA. The full line shows a fit of Eq. (§8|) with v = 0.18. The best fit 
is obtained with v — (dashed line). The dashed dotted line demonstrates the discrepancy 
between the experimental data and Eq. (|28f) when the parameter set obtained from the 3PE 
data (Fig. 3) is used. This may indicate that there are extra dephasing processes in these 
systems other than those described by the present model. 

The parameters for the heating process are summarized in Table II. There is close agree- 
ment between the estimations for 7\ presented in Sec. II and the fit values. For PMMA, 
the fit value of the diffusion time r a is about one order of magnitude larger than the estima- 
tions in Table I. For pulse intensities less than 100 nJ in Meyers' and Wiersma's experiment 
the agreement between estimated (Eq. |2[) and fitted value for T\ is less close, although the 
discrepancies are insignificant given the uncertainties in the experimental parameters for 
EtOD. 

It must be stressed that a fit to the 3PE data is also possible, if one uses the standard 
model with v = 0. However, the heating temperature we had to assume in this case were too 
high. Since v ^ introduces already a weak bending of the decay rates, it also brings down 
the value for 7\ to realistic values. We consider this observation as an indication that it is 
necessary to change the standard model in the short time regime. For longer waiting times, 
~ 10~ 3 s, the experiment explores that part of the distribution function, which becomes 
more and more hyperbolic. Hence, the bending of the decay rates should become weaker 
and, eventually, merge into the linear growth with log(t) usually observed in HB experiments. 
We believe that the deviation of the data from our theory in the millisecond regime indicates 
this transition. 
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VI. CONCLUSIONS 



We have extended Silbey's modified standard model for dephasing and spectral diffusion 
in optical experiments to the situation of a nonequilibrated phonon bath. In particular, we 
have studied the situation in which the focal volume of a laser beam in an amorphous host 
is heated to some temperature Ti, and cools down to the temperature of the helium bath 
by heat diffusion. The resulting nonlogarithmic line broadening of chromophore molecules 
embedded in the glass due to spectral diffusion has been calculated. Our motivation has 
been to investigate whether this picture explains recently observed plateaus and bumps in 
3PE data. We conclude that laser heating effects are important for 3PE experiments in the 
kelvin regime. This together with a modified, slightly bended equilibrium line broadening 
quantitatively accounts for the experimental data. However, since no studies of the 3PE 
decay on laser fluence are available up to now, the experimental evidence in favor of heating 
artifacts is indirect. A study of the fluence dependence remains an essential experiment to 
test the significance of heating effects. 
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FIGURE CAPTIONS 



FIG. 1: Experimental setup for stimulated photon echo measurements. "50/50" marks 
50 % beam splitters, DL1-3 are optical delay lines, D photodiodes, and PMT is a 
photomultiplier tube. 

FIG. 2: (a) 3PE decay rates, [l/VT^f] (t) , as a function of the waiting time, t, for ZnP in 
EtOD, from Ref. |I5|. The lines are fits to the data: ( — ) heating effects included [Eq. 
(Pop], (— • — ) no heating effects included [Eq. (]T9|)] . Effective relaxation rates used 
are R e{i = 1.2 x 10V 1 (o), R cS = 2.8 x 10V 1 (*), R e{i = 3.5 x lOV 1 (x). The 
parameters for the heating process are r a = 15 /is, Ti = 2.65 K (o), T\ = 2.8 K (*), 
and Ti = 3.3 K (x). The values of the 2PE decay rates are also given in the figure 
(+) along with the prediction of the model (dotted lines) with K v = 26.9 MHz/K 1+I/ . 
The contribution of r PLM is neglected in the fits at all three temperatures, (b) Same 
as (a) but only the data at To = 1.75 K are shown. 

FIG. 3: (a) 3PE decay rates, [1 /vrT 2 *] (t) , as a function of the waiting time, t, for ZnTPP in 
PMMA, from Ref. |H| The lines are fits to the data: ( — ) heating effects included [Eq. 
(ffiD], (— ■ — ) no heating effects included [Eq. QT9])]. Effective relaxation rates used are 
R cK = lO 1 ^- 1 (o), R eS = 1.3 x lO 1 ^- 1 (*), R eS = 2.6 x lO 1 ^" 1 (x). The parameters 
for the heating process are r a = 300 /xs, Ti = 1.5 K (o), Ti = 1.7 K (*), and T\ = T 
(x). The values of the 2PE decay rates are also given in the figure (+) along with the 
prediction of the model (dotted lines) with K v = 15.2 MHz/K 1+!y , b = 0.9 GHz, and 
AE = 9.7 cm -1 , (b) Same as (a) but only the data at T = 0.75 K are shown. 

FIG. 4: 2PE decay rate, [l/7rT 2 *] 2 p E , for ZnTPP (o) and TBT (*) in PMMA. The solid line is 
a fit of Eq. (H) to the data for v = 0.18, K v = 24 MHz/K 1+i/ , b = 15.8 GHz, and 
AE = 9.7 cm" 1 . The dashed line is a fit of Eq. @ with v = 0, K = 25.7 MHz/K, 
b = 12.7 GHz, and AE = 8.8 cm" 1 . The dashed dotted line shows Eq. with 
parameter values obtained from the fit to the 3PE data [Fig. 2(a)]: v = 0.18, K v = 
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15.2 MHz/K 1+1/ , b = 0.9 GHz, and AE = 9.7 cm" 1 . 

TABLE CAPTIONS 

Table I: Diffusion times in PMMA with p = 1.18 g/cm 3 , c = 4.6T/iJ/gK 2 + 29.2 T 3 /iJ/gK 4 
(Ref. 0), and k = (15, 20, 30, 42, 50) x 10~ 3 W/mK (Ref. |22|) at different temperatures 
T . 

Table II: Heating parameters. 

TABLES 

Table I: 



To [K] 


T a [H 


t l [ms] 


T r [ms] 


0.75 


3.10 


0.15 


2.21 


1 


4.99 


0.25 


1.99 


1.5 


10.37 


0.53 


1.84 


2 


17.05 


0.85 


1.79 


3 


43.03 


2.15 


1.75 



Table II: 





T [K] 


T est [ K ] 


If [K] r a H 


PMMA 


0.75 


1.3 


1.5 300 




1.5 


1.7 


1.7 




3 


3 


3 


EtOD 


1.75 


2.55 


2.65 15 




2.4 


2.8 


2.8 




3 


3.3 


3.3 
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FIG. 2. 
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FIG. 3. 



ZnTPP in PMMA at 0.75 (o), 1 .5 (*) and 3 K (x) 
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ZnTPP in PMMA at 0.75 K 
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FIG. 4. 



ZnTPP (o) and TBT (*) in PMMA 
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